Mon. Not. R. Astron. Soc. OOP. [TlfSTI C2005') Printed 1 February 2008 



(MN I*TeX style file v2.2) 



Magnetic Fields in Massive Star Forming Regions 



1^ 
o 
o 

Oh' 
D 



o 



> 
in 

(N 
O 

o 
l> 
o 



R. L. Curran^'^*, and A. Chrysostomou^'^ 

^School of Cosmic Physics, Dublin Institute for Advanced Studies, 5 Merrion Square, Dublin 2, Ireland 

^School of Physics, Astronomy & Maths, University of Hertfordshire, College Lane, Hatfield, Herts., ALIO 9AB, UK 

^ Joint Astronomy Centre, 660 N. A'ohoku Place, University Park, Hilo, Hawaii 96720, U.S.A. 



Accepted 



Received 



ABSTRACT 

We present the largest sample of high-mass star-forming regions observed using sub- 
millimetre imaging polarimetry. The data were taken using SCUBA in conjunction 
with the polarimeter on the JCMT in Hawaii. In total, 16 star forming regions were 
observed, although some of these contain multiple cores. The polarimetry implies a va- 
riety of magnetic field morphologies, with some very ordered fields. We see a decrease 
in polarisation percentage for 7 of the cores. The magnetic field strengths estimated 
for 14 of the cores, using the corrected CF method, range from <0.1 mG to almost 6 
mG. These magnetic fields are weaker on these large scales when compared to previous 
Zeeman measurements from maser emission, implying the role of the magnetic field 
in star formation increases in importance on smaller scales. Analysis of the alignment 
of the mean field direction and the outflow directions reveal no relation for the whole 
sample, although direct comparison of the polarimetry maps suggests good alignment 
(to at least one outflow direction per source) in 7 out of the 15 sources with outflows. 
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1 INTRODUCTION 

One of the remaining problems yet to be understood within 
star formation is the collapse of a cloud into a star, and 
the relatively slow star formation rate observed. This prob- 
lem suggests there is some form of support preventing the 
clouds from collapse (at least initially) . For instance, if typ- 
ical (M 10^~* M0) molecular clouds had no support, and 
so were in free-fall collapse, with all the mass going into 
stars, cloud lifetimes would be unrealistically short, and the 
star formation rate would be up to three orders of magni- 
tude greater than observed (jMouschoviasI [l976l ) . Thermal 
pressure is weak compared to the gravitational stresses in 
the cloud, so the support may come from the magnetic field 
that permeates the gas or the pressure of turbulent eddies - 
indeed it is likely that these two mechanisms are coupled. 

Spinning elongated dust grains in such star forming 
regions can become aligned to the local magnetic field, 
such that the grains' semi-major axes are perpendicular 
(on average) to the magnetic field lines. The net ther- 
mal emission from such grains is polarised (perpendicu- 
lar to the direction of the magnetic field lines), and can 
be used to trace the morphology of the magnetic field as 
seen projected onto the plane of the sky. Further informa- 
tion can be gained on the magnetic field via Zeeman split- 



E-mail: rcurran@cp.dias.ie 



ting fearma. Troland fc RomnevI l200lf ). which yields the 
line of sight magnetic field strength. Ion/neutral linewidths 
( Houde e t alli200d ). can also be used to gain information on 
the orientation of the magnetic field. 

In order to understand the importance of the magnetic 
field on large scales, consistent with tracing the magnetic 
field in the envelope of these cores, we present SCUBA 
850/im imaging polarimetry of 20 star forming cores. We also 
introduce a new, improved (more representative) method of 
data reduction for SCUBA polarimetry, which presents the 
errors in polarisation in a truer light. We discuss the im- 
plied magnetic field morphology across the cores, as well 
as the estimate d field strength for some of the regions us- 
ing the modified [Chandrasekha7^^^rm| ([l95^, henceforth 
CF) method. We use this estimate to indicate the impor- 
tance of the magnetic field on these scales on the support 
of these clouds, and compare the mean magnetic field direc- 
tion to the outflow directions, to see if they correspond in 
the plane of the sky. 



2 OBSERVATIONS & DATA REDUCTION 

The sample of high-mass star-formi ng regions were selected 
from the Arcetri w ater maser atlas (jComoretto et aLlll990l : 
iBrand et al] |2001^ based on the following criteria: i) well 
known and studied star-forming regions, ii) varied morpholo- 



2 R. L. Curran et al. 



gies, iii) submillimetre bright jjenness. Scott fc PadmanI 
I1995I ): and iv) observable from the JCMT in Hawaii (i.e. 
-40° ^ (5 s; +70°). 

The observations took place over several nights between 
1998 May 16 and 2000 October 10 at the James Clerk 
Maxwell Telescope (JCMT) in Hawaii. The Submillimetre 
Common User Bolometer Array (SCUBA) (jHolland et al.l 
I1999I I was used in conjunction with the polarimeter - con- 
sisting of a rotating qu artz half-wave r etard er ahead of a 
fixed wire-grid analyser; ICreaves" et al.1 (|2003l ) - which was 
mounted on the entrance window to SCUBA. The Jiggle- 
mapping mode of observation was used, involving 'jiggling' 
the secondary mirror to fully sample the beam. Sixteen dif- 
ferent positions (1 second exposures at each) are required 
to fully sample the map - these positions have a separation 
of 6.18" for the Long- Wavelength (850pim) array. The sec- 
ondary mirror performed the usual chop during the jiggle 
pattern to provide atmospheric cancellation. The nod was 
carried out over periods of around 10-20 seconds to elimi- 
nate slowly varying sky gradients. 

The polarimeter complicates the normal observing pro- 
cedure such that complete 16-point jiggle maps are required 
at specific positions of the half-wave retarder, each separated 
by 22.5°. Therefore 16 jiggle maps are observed to complete 
one cycle of the retarder. The direction of the chop-throw 
was decided based on the morphology of the target. 

The data reduction was carried out using the rou- 
tines from the SCUBA U ser Reduction Facility (SURF) 
ijjenness fc LightfootI Il998l ) to reduce the SCUBA im ages 
and routines from POLPACK (|Berrv fc Gledhillll200ll ) were 
used to reduce the polarimetry. The nod and chop of the 
telescope were corrected for and the flat-field applied to the 
observations in the standard manner. The atmospheric ex- 
tinction was calculated based on the start and end times 
of each observation: throughout the night the Caltech Sub- 
millimetre Observatory (CSO) phase monitor measures the 
T225GHZ (or Tcso), sxid a polynomial was fitted to the mea- 
sured points. This polynomial was then applied to the start 
and end times of the observation to calculate the tcso for 
the observation and converted into Tgso/jm using: 

rgsoMm = 3.99 X (tcso ~ 0.004) (1) 

for the data taken before October 2000, and: 

rgsoMm = 4.02 X (tcso - 0.001) (2) 

for the data taken during and after October 2000. The ex- 
tinction is assumed to vary linearly throughout the obser- 
vation. The airmass at which each bolometer measurement 
was made is calculated then multiplied by the zenith sky ex- 
tinction. Each data point was then multiplied by the expo- 
nential of the optical depth to give the value that would have 
been measured in the absence of the atmosphere. Bolome- 
ters which were deemed excessively noisy were switched off 
at this point in the data reduction. The sky noise was re- 
moved using bolometers that had no significant fiux from 
the source. The average (mean) fiux from these bolometers 
was assumed to come from sky emission, and was subtracted 
from all of the bolometers in order to remove the sky sig- 
nal. The instrumental polarisation was removed in the usual 
manner. The data were re-gridded using a Gaussian weight- 
ing function, with the scale set to 7" (half the beamsize). 



Table 1. Tabic of flux calibrators and flux calibration factors. 



Date 


Flux Calibrator 


FCF 






Jy/arcsec^ /V 


19980516 


Uranus 


2.35 


19990620 


Uranus 


2.02 


19990705 


Uranus 


2.21 


19990706 


Saturn 


1.81 


20001006 


Uranus 


1.96 




CRL2688 


1.72 


20001010 


CRL2688 


1.71 



The pixel size was set to 6.18" , matching the jiggle pattern, 
in order to calculate the polarimetry vectors accurately. 

POLPACK packages are then used to calculate 
the Stokes paramete rs by fitting the following curve 
jSparks fc Axonll 19991) to the data: 

l'k^^{I + e{Q cos 2cj,k + Usm2cj,k)) (3) 

where 7^ is the expected intensity in image k, t is the wire- 
grid analyser transmission factor, e is the analyser polarising 
efficiency factor and 0^ is the effective retarder position an- 
gle after correction for the parallactic angle for image k. 
The polarisation percentage and position angles of the vec- 
tors are then calculated from the Stokes parameters, without 
any binning. The catalogue was clipped such that noisy po- 
larisation vectors were not included. The clipping used was 
7 > and dP < 0.75. For 3% polarisation, this represents 
a maximum position angle error of 7.15°. It was chosen to 
clip the vectors on polarisation errors instead of signal-to- 
noise, as clipping on the latter would result in disposing of 
points where the polarisation is low or zero, both of which 

are perfectly valid measurem ents. 

Finally, flux calibration (|jenness et al.l [20021 ) was car- 
ried out using observations of Uranus, Saturn and CRL2688. 
The fiux calibration factors (FCFs) are listed in table [T] 
These FCFs were applied to the data, with any objects ob- 
served over more than one night having the data for each 
night reduced separately, then the flux calibrated images 
co-added to create the final image. 

2.1 Pixel Size 

Since the introduction of SCUBA on the JCMT, it has be- 
come "standard" to reduce 850pim data with a pixel spacing 
of 3.09" . This is mainly due to a 3.09" jiggle-step being re- 
quired to fully sample the beam at 450/im. This allowed for 
easier comparisons between 450/im and 850^m data. The 
method of data reduction naturally continued when the po- 
larimeter was introduced. 

Careful inspection of data reduced by this method (pixel 
size set to 3.09"and with 2x2 binning of vectors) and re- 
duced using a pixel size of 6.18"with no binning, revealed 
significant differences in the number of vectors. The polari- 
sation errors in the latter reduction were higher than those 
in the first reduction method, resulting in a decrease in the 
number of vectors selected after error clipping. 

The change in the errors in polarisation are due to the 
binning. When the pixel size is set to 3.09", four pixels are 
produced for every one when the size is 6.18". A vector is 
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then produced for each pixel, with the vectors for the smaUer 
pixels being binned together — adding the polarisations to- 
gether in quadrature to calculate the error. This erroneously 
reduces the resultant errors. 

The observing mode used means that one independent 
polarisation measurement is taken at steps of every 6.18" , 
the results of which can be superimposed over an intensity 
image created with 3.09" pixels. 



3 CORE PROPERTIES 
3.1 Core Masses 

The mass of the star-forming cores can be calculated by 
assuming: 



de = 



28.6° 



(6) 



MTota 



(4) 



where g is the gas-to-dust ratio, Sv is the flux, d the distance, 
K,u the absorption coefficient at frequency u, and -B,y(Tdust) 
is the Planck function for frequency v at a. temperature of 

Tdust. 

Using a gas-to-dust ratio of 100: 1 l|Hildebrandlll983h . 
and an absorption coefficient of 0.15 m^ kg~^ estimated from 
lOssenkopf fc Henningj (|l994l ') based on a number density of 
ilh = 10® cm~^, thick ice mantles and a formation timescale 
of 10® years, mass estimates for the cores were calculated 
(listed in table [2j along with the other parameters derived 
from the observations). 

The mass of these objects are subject to uncertainties 
within the used parameters, especially any errors in mea- 
suring the distance to t he core. The gas-to-dust ratio, which 
may be as low as 45:1 (|McCutcheon et al.lll995l ) also intro- 
duces an error that may be up to a factor of 2. The ab- 
sorption coefficient at 850Mm still has not be en determined 
precisely l|Hildebrandll 19831 : IChini et al.ll2001^■ although the 



value adopted in this paper from lOssenkopf &: Henning 



1994) agr ees with the values determined bv lBianchi et al 



20031 ) and|v isser et a.l(|2002l ). The masses we determine in 



this paper are therefore the upper limits for the temperature 
of the cores. 



3.2 The Magnetic Field Strength 

The CF relation can be used to obtain the plane of the sky 
average magnetic field strength. This method is based on 
equipartition and the ability of the magnetic field to retain 
straight field lines under the infiuence of turbulence. The 
plane of the sky average field strength can be calculated 



G, (5) 

where p is the mean density (g cm~^), (Tt,,^^ the line of 
sight velocity dispersion (cm s~^), ag is the dispersion in 
polarisation position angles (measured east of north) and is 
corrected for measurement errors {ag = crj^;,jjg^j;,d — (Terror) 
where 6 is in radians, and / is a correction factor found to 
be ~0.5 (|Heitsch et al.ll200lD . 

Errors on the position angle of the vectors are calculated 
based on the polarimetry signal-to-noise: 



where d9 is the error in the position angle and Sp is 
the signal-to-noise in the polarisation. Therefore in regions 
where the polarisation percentage is low (e.g. across the 
main core of DR21(0H)), then ameas ~ cTen- resulting 
in misleadingly small dispersion angles, and very strong 
magnetic fields. By only selecting the vectors which have 
a signal-to-noise (in polarisation) of > 3, this problem is 
avoided, therefore all of the magnetic field strength esti- 
mates in this paper are based on vectors selected in this 
way. Please note, however, as mentioned earlier, the vectors 
shown in the figures are not selected in this way, as in terms 
of magnetic field morphology, polarisation nulls and low lev- 
els of polarisation are valid detections. 

The CF method should be used with caution when cal- 
culating the magnetic field strength in the plane of the sky. 
The 15" beam-size of the JCMT means that the small- 
scale tangling of the magnetic field field can occur within 
the beam, and so the measured vectors only represent the 
net magnetic field direction, leading to over-es timates of the 
field strength. Modelling studies of this effect l|Heitsch et al.l 
120011 ) have lead to the introduction of a correction factor /, 
which has been found to be ~ 0.5 (see eq. [5]). Also, any un- 
derlying magnetic field morphology (e.g. intrinsic field cur- 
vature) has not been accounted for in calculating the dis- 
persion in position angles of the polarimetry vectors (i.e. 
this technique assumes a priori that the magnetic field is 
uniform) . 

Uncertainties in our estimations of the magnetic field 
strength arise from calculating the density of the cores, 
which incorporates the errors involved in calculating the 
mass, they therefore represent the upper limits of the mag- 
netic field strength. There are also errors in calculating the 
volume of the core which contribute as a spherical geom- 
etry has been assumed for each core (except where other- 
wise stated), although with no density tracer information, 
we lack the ability to modal the three-dimensional structure 
of the cores, therefore the upper limits stated are solely for 
spherical geometry. The velocity of the gas within the core 
introduces another error as a FWHM of ~ 2 kms^^ has 
been used but it may be anywhere between 1 kms^^ and 
3 kms"^ l|Brand et all [2001: Thomp son et al. 2004). Mea- 
surement errors are also introduced by the angle 6, although 
these are relatively small in comparison to the other errors 
stated. 



4 INDIVIDUAL CORES 
4.1 Cepheus A 

Cepheus A is a well-known st ar forming region. Located 
at a distance of ~730 pc (|Blaauw. Hiltner fc JohnsonI 
Il959ll . it is the closest star-forming region in the sample. 
The region contains many signatures of massive star 
formation such as: a) a sub-mm source (observe d here) 
with a total luminosity of 2.5 x 10* Lq (Koppenaal et al.l 
1 19791 : lEvans et al.lll98l|). b) a cluster of 14 compact radio 
continuum sources (iBeichman. Becklin fc Wvnn- Williams! 
1 19791 : iHughes fc Waterloutlll984l). c) clusters of H2O and 
OH masers (|Lada et all Il98ll : iMigenes. Cohen fc Brebnerl 
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04 02 20:39:00 38:58 



Right ascension 

Figure 1. Final result of the data reduction method outlined in section 2. Epoch J2000. The greyscale is the total thermal continuum 
emission at 850/ini, with the B-vectors overlaid to represent the plane of the sky magnetic field direction. The polarisation percentage 
scale is presented in the lower left of each image, with the total intensity scale to the right of each image. 
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Table 2. The sources with their observed parameters. Magnetic field strengths are not calculated for all cores, e.g. DR21(OH) Main. 



The vectors across the 


cores which have no field 


strength estimates have low polarisation signal-to- 


noise 


(see section 3.2). 
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DEC 


Distance 


Temp 


Flux (850Atm) 


Mass 




n 
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Temp. Ref." 




(hh;mm:ss) 


r " ) 


(kpc) 


K 


(Jy) 


(Mq) 


cm 


-3 


mG 




Cepheus A 
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9 



° Where no reference is given, temperatures are assumed to be 20 K, consistent with the coldest measured temperature in this sample, 
DR21(OH). This therefore leads to upper limits for the masses of these objects. 
For a cylinder of radius 0.33pc (22.7" ) and height 0.93 pc (64" ). 

For a 40" aperture, instead of a 60" , to stop contami natio n be tween sources. 

Temperature references a re: (1) BotineUi & WiUiams (20^); (2 )lManeum. Wootten fc Mundvl lll99^: f3}lvan der Tak et al.l ||2000^ : (4) 

bl4200i^ 



McCutcheon et al.l lll995h : (5) ISandcll &c Weintraub 
Harvey et al.l l ll977l) 



, - . . --T 

(6) iThronson et al.l l|l980l '): (71 ISandelll l|l999l ) : CS) IWu et al.l l|2005^ : (9) 



Il992l'). d) exte n ded s h ock-excited H2 em i ssion at 2.12/im 
ijBallv fc Lan3 llQSl iDovon fc Nadeaul 1 19881 ) and e) 
highly blueshifted Herbig-Haro objects with large 
proper motions pointin g away from the actiyity centre 
dLenzen. Hodapp fc Solj 119841: iHartigan fc Ladal 1 19851 : 



iHartigan et al 



a,pp &c 501 

.1 1 19861 : iLenzenll 19881 ). The well known outflow 
in this region has a complex morphology. On scales of 
~ 1 pc, the outflow is in the east-west direction, but 
on smaller scales, th ere is an ionised jet t hat is oriented 
northeast-southwest. [Narayanan fc Walked (|l996) argue 
that this may be due to multiple outflow episodes. They 
hypothesise that the older (3-20 xlO*yr) high- velocity 
outflow is oriented along the east-west direction, while the 
extremely high velocity lobes of the current flow (^ 5 x 10^ 
yr) extend to the northeast and southwest. On even smaller 
scales, there is evidence for a thermal radio jet with position 
angle ~ 48° , which is perpendicular to a band of water 
mase r positions that ar e thought to trace a circumstellar 
disc l|G6mez et al.lll999l ). 

The magnetic fleld in Cepheus A has been stud- 
ied recently using Zeeman effect of H2O masers 
llVlemmings et al.ll2006l) and OH masers (IBartkiewicz et al.l 
I2OO5I ). Vlemmings et all (|2006l ) find that the maser emission 
associated with the HW2 high-mass young stellar object 
in this region, which is argued to trace the circumstellar 
disc, has a field strength of 100-600 mG. They find other 
meisers further from HW2 to have magnetic field strengths 
of 30-100 mG. In all cases the magnetic field pressure is 
calculated to be similar to the dynamic pressure, indicating 



the magnetic field is strong enough t o cont rol the outfiow 
dynamics of HW2. IBartkiewicz et al.l l|2005l ) find that the 
magnetic field is extremely ordered on arcsecond scales, 
pointed away from us in the east, and toward us in the 
west. Their linear polarisation measurements from these 
masers reveal that the direction of the magnetic field in the 
plane of the sky is ~ 53°, or roughly parallel to the radio 
jet axis. 

The SCUBA polarimetry of Cepheus A (Fig[T^) reveal 
magnetic field vectors that are ordered around a position 
angle of 50° in the northeast of the observed region, and 
ordered ~ 40° in the southwest of the region. There is a 
band of depolarisation across the core which has a position 
angle of ~ 135° perpendicular to the inferred magnetic fleld 
direction. This plane of the sky magnetic field morphology is 
roughly parallel to the outfiow/radio jet axis that has been 
observed on similar/smaller scales. This pattern of magnetic 
field vectors may be interpreted as a 'pinched-in' or hour- 
glass morphology, with the field twisted towards the centre 
of the source, where the depolarisation occurs. 

CF calculations for the plane of the sky magnetic 
field strength lead to estimates of ~ 6 mG. This is the 
strongest magnetic field calculated for the cores in this sam- 
ple, but is smaller t han the field stren gths calculated from 
the maser emission (| Vlemmings et al.l l2006). This suggests 
that the magnetic field is weaker on larger scales, becom- 
ing stronger on smaller scales, closer in to the protostar. 
The data presente d here , in c onjunction with the findi ngs of 
IBartkiewicz et ahl (|2005l ) and I Vlemmings et al.l l|2006l ). lead 
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to a magnetic field that, in the plane of the sky, is northeast- 
southwest in direction, with the magnetic field in the north- 
east pointing away from us, and in the southwest pointing 
toward us. 

4.2 DR21(OH) 

DR21(0H) is part of the star forming complex W75, lo- 
cated in the Cygnus X region, at a distance of ~ 3 kpc 
(jCampbell et al.lll982l ). Its relative close proximity has aided 
in the fact that it is one of the most studied star form- 
ing regions in the Galaxy. The DR21 H ll region is ~ 
3' to the south of DR21(0H) and contains a cluster of 
late type O stars, and one of the most intense outflows 
known. DR21(0H) - also known as W75S or W75S(OH) 
- unlike DR21, consists of young stars in the process of 
forming, and as such offers a glimpse of an earlier stage 
of evolution of star formation. Previous continuum stud- 
ies have shown that DR21(OH) is made up of three com- 
pact continuum sources, DR21(0H) Main, DR21(0H)S and 
DR21( 0H)W, all of which seem to be actively forming 
stars ijMangum. Wootten fc Mundvl Il99ll ). Higher resolu- 
tion studies confirmed that DR21(OH)Main i s in fact com- 
posed of two smaller cores, MMl and MM2 l|Woodv et al.l 
1 19891 ). 

The magnetic field of DR21(0H) is reasonably well 
studied, and there are Zeeman estimates of the line of sight 
magn e tic field strength for MMl and MM2 (jCrutcher et all 
ll999l '). lLai et all l|2003bl l used BIMA to gain information on 
the magnetic field via polarimetry of both dust and CO, in 
order to map the magnetic field morphology in the plane 
of the sky. Comparisons of ion/neutral linewidths have also 
been used to est ablish a three-dim ensional impression of the 
magnetic field (|Lai et al.l l2003al ) . The morphology of the 
field at the resolution sampled by BIMA indicates that the 
magnetic field is ordered, thus implying a strong field. The 
magnetic field is estimated to have a strength of ~ 0.4 mC 
with an inclination of 36° to the line of sight, and a pos i- 
tion angle of 105° in the plane of the sky ()Lai et al.ll2003al ). 
The BIMA data only measure the polarised dust emission 
in patches, and do not reveal the magnetic field morphology 
throughout the whole region. 

The data presented here (Fig[l]3) include our observed 
data as well as some archival d ata for DR21(QH ), which were 
observed on 2002 October 2 ijVallee fc FieedEoOS ). which 
effectively doubles the time on source for this target. The 
archival data were subjected to our method of data reduc- 
tion as detailed earlier. The data reveal that the DR21(0H) 
region is composed (at this resolution) of one main core 
DR21(0H) Main. There are also two fainter cores close to 
the main core, one to the s outh-west, DR21(OH)W and 
one d ue south, DR21(OH)S. iMangum. Wootten fc Mundvl 
l|l992h identified another core, DR21(0H)N, in NH3 emis- 
sion, although in the SCUBA data presented here it appears 
to have the morphology of a ridge extending northwards 
from the main core. There are also molecular outflows asso- 
ciate d with the main core emanating in an east- west direc- 
tion (|Lai et al.ll2003bl ). 

Our calc ulated mass of the Main core is somew hat larger 
than those of lMangum. Wootten fc Mundvl (119911 ). but they 
calculated the individual masses of MMl and MM2 using 
OVRO interferometer measurements of the dust emission. It 



is to be expected that our calculations reveal higher masses, 
given that they are based on single dish measurements and 
so include more diffuse dust on larger scales. Comparisons of 
our derived total masses with those of,yallee fc Fiegc (200^ 
reveal our masses are much higher than their estimates, how- 
ever they have assumed a temperature of 100 K for both the 
Main and northern source. We also use larger apertures for 
our calculations. Once these two factors have been taken 
into account, the mass estimates are consistent. 

The polarimetry indicates that the magnetic field is 
ordered across the ridge, parallel to the outfiow axis. The 
percentage polarisation also remains stable (at ~ 3%) 
across the ridge. Across the main core the percentage po- 
larisation drops, most apparently to the northeast of the 
core, coincident with MMl. T his is observed in numer- 
ous other cores (for example, 'Chrvsosto mou et al.l I2OO2I : 
Matthews & Wilson 2002; Davis ct al. 200_^) and could be 
due to the magnetic field twisting within the JCMT beam, 
the grains becoming more spherical in regions of high den- 
sity, or the grains being less efficiently aligned in regions of 
high density. 

To the south of the main core, the vectors across 
DR21(0H)S are more dispersed in position angle, such that 
to the southwest of DR21(OH)S, the vectors have a position 
angle of ~ 135°, changing to ~ 90° northeast of the source. 
The vectors across the south-western core are ~ 90°, the 
same as across the ridge. 

The polarime t ry dat a, in general agree with the findings 
of lVallee fc Fiegel l|2006l '). however we see a smaller disper- 
sion in position angles throughout the region, which may be 
due to a combination of higher signal-to-noise observations 
and a more careful data reduction. The polarimetry agree 
with the findings of iLai et al.l l|2003bl ). with the polarisation 
nulls coincident with MMl, where Lai et al. found very lit- 
tle polarisation from dust. The overall change in direction of 
the magnetic field across DR21(0H)Main also agrees with 
the BIMA data - both the dust and the CO polarimetry 
- indicating that the magnetic field stays ordered on both 
large and small scales. The observed polarisation nulls across 
MMl may indicate that either MMl has a twisted magnetic 
field in comparison to MM2, or that MMl is more centrally 
condensed than MM2. Both of these would be consistent 
with MMl being the more evolved of the two cores. 

The polarimetry across the ridge (northern source) of 
the DR21(OH) region indicate an ordered field, which in it- 
self implies a strong field. The field is (in the plane of the 
sky) perpendicular to the north-south ridge, yielding the 
possibility that collapse has occurred along the field lines. 
CF calculations reveal the magnetic field strength is of the 
order ~ 1 mC in the plane of the sky. This is compara- 
ble to the field strengths previ ously calculated for th e two 
sources in DR21 (OH)Main by ICrutcher et al. I l|l999l ) and 
iLai et al.1 (|2003bh . which were ~ 1 mG in the plane of the 
sky from BIMA observations and ~ 0.5 mC in the line of 
sight from Zeeman measurements. This indicates that the 
magnetic field is not only morphologically uniform through- 
out this cloud, but the strength is also re asonably uniform 
on different scales in the cloud. iVallee fc Fiege (2006) esti- 
mated the strength of the magnetic field across the region 
as a whole (finding 780/xG), across the main source (finding 
780/iG again) and across the ridge to the north of the main 
core (~ 200/iG) . Our estimate of the magnetic field strength 
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across the ridge is consistent with their estimate across the 
entire region, but much greater than their estimate across 
the ridge, although they only use the four closest vectors to 
the source peak to calculate their field strengths, whereas 
we use all of the vectors which fall within the aperture used 
for the mass and density calculations. 

4.3 GGD27 

GGD27 is located within Sagittarius, on the southwestern 
edge of a dark lane that runs in a northwest-southeast 
direction. It's kinematic dista nce has been calculated to 
be 1.7 kpc (|Rodrfguez et al' via the velocity of 

observed CO. It is a well known site of star forma- 
tion, and located within the SCUBA submillimetre core 
(marginally reso lved by the JCMT b eam) are several in- 
frared sources (jStecklum et al.l 1 19971 ). There is a large- 
scale CO north-south outflow (bl ue lobe in the north ) 
from an embedded source (IRS2) lYamashita et al.1 Il987l ) 
as well as a CS disc elongated in t he east-west direction. 
iMartf. Rodrfguez fc ReipurthI (|l993l ) discovered one of the 
largest radio jets emanating from this region - it's source be- 
ing coincident with IRAS18162-2048, an extremely luminous 
star (2x IO^Lq). They found the radio jet has a position an- 
gle of 21° , extends 5.3 pc, and is extremely collimated, with 
an opening angle of just 1°. 

The polarimetry data of GGD27 (Fig[T};) imply a mag- 
netic field that has a north-south direction, roughly paral- 
lel to the molecular outfiow. To the north of the core, the 
polarimetry show that the magnetic field 'fans out' - the 
magnetic field vectors to the north-east of the core have a 
north-east orientation, whilst the vectors to the north-west 
of the core have a north-west orientation. This could be an 
indication that the field is pinched in closer to the core, 
although no depolarisation is observed across the core, sug- 
gesting that the magnetic field is not sufficiently twisted. 
The plane of sky magnetic field strength is estimated to be 
~ 0.2 mG. 



4.4 GL2136 

Near-infrared studies of GL213 6 (|Minchin et al.1 Il99ll : 
iKeistner. Weintraub fc AsTainlll993 ) indicate that there is a 
circumstellar disc or torus roughly at a position angle of ^ 
45° . A kinematic distance of 2 kpc has been calculated by 
iMenten fc van der TakI l|2004h . Both water and OH maser 
emission has been observed, with analysis of the left and 
right circularly polarised components of the 1665 GHz fea- 
ture leading to line of sight ma gnetic field estimates of 1 mG 
ijMenten fc van der TaM l2004h . CO observations have re- 
vealed a massive (50Mq) molecular outflo w from the sou rce, 
perpendicular to the observed disc (Kast ner et al.|[l9 94i). 

The submillimetre continuum image of GL2136 (Fig[ljl) 
reveals a marginally resolved object with a small extension 
to the southwest. The magnetic field vectors are mainly in 
the east-west direction, however in the south of the source 
they are more northwest-southeast orientated. The average 
position angle is 102° , around 30° away from the outfiow 
direction. There is no depolarisation towards the centre of 
the core. CF calculations reveal plane of sky magnetic field 
strengths of ~ 0.3 mG. This is less than half the value for 



the line of sight component measured from maser emission, 
therefore, either the magnetic field is weaker on large scales, 
or, if we envisage the magnetic field wrapping around the 
outflow axis (which is close to the plane of the sky - see ta- 
ble |3}, it must have a large toroidal component to produce 
the large line of sight component of the field strength - in- 
deed, this may explain the relatively large (^30° ) difference 
between the outfiow axis and the magnetic field direction. 

4.5 GL2591 

GL2591 is a well-studied source, although its distance is 
still quite uncertain. Distance estimate s range from 1 kpc 
(iMozurkcwich. Schwartz fc SmithI 'l986) to greater than 2 
kpc ((Merrill fc Soifer 19?! ). The majority of estim ates are 
between 1 and 2 kpc (e.g. IWendker fc Baar3ll974l ). there- 
fore a mid-range di s tance of 1.5 kpc is adopted here. 
iHasegawa fc Mitchelll l| 19951 ) studied the molecular outfiow 
in detail, finding a small-scale, well-coUimated fast outflow in 
an east-west direction superposed on the previously known 
northeast-southwest, large-scale, less-coUimated slow out- 
fiow; 

iHutawarakorn fc CohenI (|2005h studied the magnetic 
field using Zeeman splitting of OH masers. The masers 
formed an elliptical shape perpendicular to the outflow axis, 
which was suggested to be a molecular torus of radius ~ 750 
AU, inclined at 55° to the line of sight. The magnetic fleld 
was found to range from -1.6 to +3.8 mG, reversing direction 
on opposite sides of the disc, indicating a toroidal compo- 
nent. The linear polarisation vectors were found to be both 
parallel and perpendicular to the outfiow direction. 

The SCUBA data (Fig. [1^) show a marginally resolved 
object, and the magnetic field vectors have, in general, an 
east- west direction across the core, parallel to the small- 
scale (90" X 20" ) outfiow. The western side of the core has 
a polarisation null, and to the south of this, the mag- 
netic field appears to be in a southeast-northwest direc- 
tion. To the eastern side of the core, the vectors begin to 
curl up towards the northeast. The average position an- 
gle of the polarisation vectors is 95° aligning with the 
outflow axis. The polarisation null to the west of the ob- 
je ct suggests a twisted fie l d, ag reeing with the findings 
of IHutawarakorn fc CohenI (|2005l ). The plane of sky field 
strength is esti mated to be ~ 0.5 mG. When co mpared with 
the estimates o f IHutawarakorn fc CohenI (|2005h . which indi- 
cate a strong toroidal field in the disc, the data imply that 
the magnetic field, once again, becomes weaker on larger 
scales. 



4.6 GL437 

GL437 is located at a distance of ~ 2 kpc 

dArauifla fc GoldsmithI Il984h . It is composed of a 
compact cluster of young stars (including B stars), 
and a refiection nebulosity which is centered on WK34 
(j Weintraub fc Kastneil[l996l ). There is also a broad molec- 
ular bipolar outfiow t hat extends ~ 1 pc and is orientated 
roughly north-south (|G6mez et al.l (19921) . with the south 
lobe being blueshifted. The SCUBA data (Fig. ^) reveal 
that this source is an elongated core orientated in the 
northwest-southeast direction. The polarimetry is aligned 
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such that the magnetic field appears to be projected 
perpendicular to the ridge of gas and dust forming the 
elongated source, roughly parallel to the outflow axis. There 
is a polarisation null to the southeast of the core, and the 
polarisation percentage changes seemingly randomly across 
the core. 



4.7 IRAS 20126+4104 

IRAS 20126+4104 is located in the Cygn us-X region, at 
a distance of 1.7 kpc (IWilking et al.lll989l). It is classified 
as a high mass protostellar object (HMPO: ISridharan et al] 
|2002| ) — a lthough at a later sta ge of evolution than those dis- 
cussed in ICurran et all (|2004l ) — it is the high mass equiv- 
alent of a class O/I protostar. It has a well known molecular 
outfiow and an ionised jet. The CO outflow is roughly north- 
south in orientation (position angle of 171°) with redshifted 
gas in the south and blueshifted gas in the north. The jet, de- 
tected by emissio n knots of H2 and [Sll] has a position angle 
of 117° (.Shepherd et al.ll2000l ). [Shepherd et all (|2000l ) con- 
clude that the most likely interpretation for this is for the 
coUimated jet to be precessing thr ough an angle of ~45°. 
More recently, iLebron et al.l l|2006l ) studied the kinematics 
of this region and concluded that whilst a precessing jet is 
supported by the data, multiple flows driven by independent 
sources cannot be ruled out. 

The polarimetry data (Fig.[T^) reveal that the projected 
magnetic fleld has, in general a northeast-southwest orien- 
tation. This is roughly at an angle of ~ 40° to the ionised 
jet, and almost perpendicular to the CO outflow, which may 
be evidence of a helical field. There are several polarisation 
nulls to the north of the core, away from the intensity peak. 
The plane of sky component of the magnetic field is small - 
only 0.1 mC 



4.8 IRAS 20188+3928 

IRAS 20188+3928 is located in the C ygnus region at a n 
uncertain distance between 0.4-4 kpc l|Little et aLlllQS"^ . 
All further calculations in this paper assume a maximum 
dist ance of 4 kp c , and so all values quoted are upper lim- 
its. iLittle et al.l (|l988l ) observed this region in CO and 
HCO^, and found a CO bipolar outflow, which has a north- 
southwest direction. The redshifted gas is to the north, and 
the blueshifted gas is in the southwest. The SCUBA data 
(Fig. [T}i) shows a curved ridge of gas and dust extending 
from the north of the core, curling around to the west. The 
polarimetry (in general) indicate that the magnetic field 
lines are east- west in orientation, projected onto the plane of 
the sky. Towards the southwest of the core, the polarimetry 
vectors begin to deviate from the east-west orientation to- 
wards a northeast- southwest orientation, aligning with the 
outfiow in the southwest. The degree of polarisation is mea- 
sured to be (on average) lower across the ridge than across 
the brighter core. A plane of sky field strength component 
of ~0.2 mC has been calculated. 



4.9 L1287 



cloud, orientated in a northeast-southwest direction. At the 
centre of the outfl ow there is a ver y cold IRAS source 
IRAS 00338+6312 (|Yang et al.l liggil V The submillimetre 
data (Fig. [TJ) reveal a marginally resolved object. There is 
little polarimetry (2 vectors) as the integration time for this 
object was only 0.85 hrs, and so the errors in polarisation 
are large. 



4.10 MonR2 



Monoceros R2 is a s ite of ongoing star format ion, with a 
compact Hll re gion llWood fc Churchwelll Il989l) and H^O 
and OH masers iDownes et al.ll 19751 : iKnapp fc Brownlll976l '). 



L128 7 is a dark cloud at a distance of 850 pc l|Yang et al.l 
I1991I ). There is an energetic outflow associated with the 



The outflow in this region is one of the largest known, 
with a total extent of 6.8 pc (assu ming a distance of 
950 pc; iRacine fc van de Berghl Il970l ). This outflow ex- 
tends to the north, east and south-west o f the region. 
ICiannaJtopoulou. Mitchell fc Hascgawa! (|l997l ) assume that 
the large scale north-southwest outflow is old, and that the 
source of this outflow is now inactive, as the lobes are not 
well coUimated. The easterly outflow is thought to be a rel- 
atively young, small-scale outflow, unresolved, from another 
source within the region. The polarimetry of this region 
(Fig-HJ) indicate that the magnetic field aligns with the out- 
fiows - the vectors are oriented north-south in the north of 
the region, east- west in the east of the region and northeast- 
southwest in the southwest of the region. The morphology of 
the magnetic field is complex, with a plane of sky magnetic 
field strength of ~ 0.2 mC 



4.11 NGC6334A 

NGC6334A is a giant Hll region/molecular cloud complex 
and is l ocated close to the galactic plane at a distance of 
1.7 kpc (|Neckellll978h . NGC6334-submm has two cores, one 
in the east, and a brighter core (probably multiple, but not 
resolved within the JCMT beam) in the west of the image. 
There is a CO bipolar outflow associate d with NGC6334A , 
which lies almost in the pla ne of the sky llSarma et al.|[2000h 
in a north-south direction l|De Free et al.l 1995f ). There is a 
ridge of material between the two cores in which the po- 
larimetry (flg. [T]<) suggests that the magnetic fleld is paral- 
lel to the ridge. The magnetic field across the faint core has 
a southeast-northwest direction, which curls round to east- 
west across the ridge. The magnetic field across the brighter 
main core is more complex, in the northeast of the core, 
the magnetic field has a northeast-southwest direction, but 
towards the southwest of the core, it is directed northwest- 
southeast, and in the very southwest of the region there is 
a polarisation null. The abrupt change in the orientation of 
the magnetic field may be due to two cores being present, 
although unresolved in the JCMT beam. Previous observa- 
tions of OH maser emission have yielded estimates of the 
line of sight magnetic field strength of ^ 0.35 mG toward 
the s ource, with maximum strengths of 0.5 mG (iMavo et al.l 
l2004h . The polarimetry presented here lead to estimates of 
the plane of the sky component of the field strength to be ~ 
0.9 mG, consistent with the previous line of sight measure- 
ments. 
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4.12 RCrA 

The molecular complex Coronae Austral is is ~ 129 pc away 
from the Sun (jMarraco fc Rvdgrenlll98lh . and is dominated 
by the centrally condensed core centred near the emission 
line star R Cr A. Nutter, Ward-Thompso n fc Andrei (120061 ) 
have recently studied this source in the submillimetre, and 
reveal three submillimetre peaks within the main source, 
named SMM-1A,B and C. SMM-IA is located in the south- 
eastern part on the core, with SMM-IB and C being located 
in the north-east and northwest of the core. 

R Cr A is one of the more evolved sources in the 
polarimetry sample, and i s classified as a Herbig Ae star 
IIMarraco fc Rvdgren|[l98ll ). There is a CO bipolar outflow 
with a position angle of ~90° in the region of R CrA 
jWalker. Lada fc Hartiganlll985l) but more recent molecu- 
lar line mapping by Anderson et al.l (|l997l ) has cast doubt 
on whether R CrA is the driving source. The polarimetry of 
the R CrA region (fig. [TJ) indicate that the magnetic field 
is roughly parall el to the direction of the outfiow (see also 
IClark et all 120001 ). There are also polarisation nulls across 
this source, in two regions — one in the north of the source, 
and one in the south. A plane of sky field component of 1 
mC has been calculated. 



4.13 S140 

S14 is located in the Ceph eus ring, at a distance of 900 
pc iPreibisch fc Smith! [200^ ) . S140 has two outflows, one 
with a position angle of ~160°, which is bipolar in nature, 
and anoth er sm aller scale outflow which has a position angle 
of ~20° (jPreibisch fc Smith 20 0^). The higher resolution 
K'-band data of IWeigeh et all l|2002l ) show arc-like struc- 
tures protruding from the northeast of the source which are 
proposed to trace outflow cavities carved out by material 
flowing away from S140 IRSl. The resolution of their data 
is 240 milli-arcsec and covers an area of 13" x 21" - ap- 
proximately one beam width of the SCUBA data presented 
here. The SCUBA data (fig. [TJn) show that the magnetic 
field vectors are ordered, in a north-south direction to the 
western side of the source, whereas in the east, the mag- 
netic field is east-we st orientated. OH Zeeman observations 
(jBaudrv et al.ll 19971 ) have revealed line of sight field esti- 
mates of -1-2.8 mC. CF estimates of the plane of sky compo- 
nent of the field lead to strengths of ~ 0.4 mC, which either 
indicates the magnetic field is mainly in the line of sight, or, 
that the field is stronger on smaller scales. There is a faint 
ridge of gas and dust extending from the east of the source, 
curling northwards, which the magnetic field vectors seem 
to follow, running parallel to it (in the plane of the sky). It 
may be possible that the arc-like structure in the SCUBA 
image is relate d in some way to th e smaller scale outflow 
cavities seen bv lWeigeh era!(|2002l ). 



dust seemingly connecting the two cores. The polarimetry 
appears to be almost randomly distributed, with several po- 
larisation nulls on the northern core, and to the north and 
south of the southern core, although there does not seem to 
be a relation between the intensity and polarisation percent- 
age. The large scatter of the polarimetry vector position an- 
gles may suggest that the magnetic fleld is weak across this 
region. It may also be explained if the magnetic field was 
predominantly in the line of sight (the outflow is mainly in 
the line of sight - see table [3]), which can cause random po- 
larisation patterns and/or low polarisation percentages. The 
plane of sky component of the field strength is calculated to 
be 0.1 mC, which is one of the weakest measured for this 
sample. 



4.15 S157 

S157 is a diffuse nebula located towards the Cassiopeia- 
Perseus arm at a distance of ~ 2.5 kpc, and is surrounded 
by Hll regions and young open clusters (iShirlev et al. 200^). 
The submillimetre data (fig. [1})) show the region is only 
slightly more extended than a point source, with bright sub- 
millimetre emission extending southwards from the main 
core. The polarimetry of the core show that the magnetic 
field vectors have an east-west direction to the south of 
the core, but in the north, the vectors are aligned roughly 
northwest-southeast. There are two regions of null polari- 
sation — one in the northwest of the source, the other in 
the southeast. The plane of sky field strength component is 
estimated to be ~ 0.2 mC. 



4.16 W49 

W49 is in the galactic plan e at a distance of 11.4 kpc 
ICwinn. Moran fc Reidlll992l ). and is one of the most lumi- 
nous Hll regions in the Galaxv lsinith et al.l (~ 10^, wf^- 
There are two bright cores in this region (fig. [T]d) - one in 
the northwest, W49N, and one in the southeast, W49SE, 
with a third fainter source along the ridge in between the 
brighter two cores, W49E. This ridge differs from those pre- 
viously mentioned, as it is much more extended and less 
concentrated (fainter). There is a CO bipolar outflow from 
W49N, with the redshifted lobe to the north and blueshifted 
to the south. The outflow is almost in the line of sight 
(iScoville et al.|[l986l ). The polarimetry vectors are aligned 
such that the magnetic field lines run from one core to the 
other, parallel to the ridge. Across the W49N core, the mag- 
netic field (plane of the sky component) is estimated to be 
less than 0.1 mC. The outfiow is in the line of sight, how- 
ever, and so this may suggest that there is a larger line of 
sight component to the magnetic field. 



4.14 S146 

S146 is located at a distance of 5.2 kpc l|Wu et al.ll2005l ). A 
bipolar molecular outflow, in a north-south direction (north 
lobe blueshifted), driven by a star of spectral type 06.5 or 
earlier (required to account for the ionisation of the Hll re- 
gion). In the submillimetre data (flg.[l]i) there are two cores 
in a north-south configuration. There is a ridge of gas and 



5 DISCUSSION 

5.1 Comparison of Outflow &: B-field directions 

Theoretically, magnetic fields play an important role in the 
launching and collimation mechanisms of outflows, and so 
it would be of interest to see if any relationship between 
the observed magnetic field direction on these large scales, 
and the jet/outflow axis exists. Previous studies of the 
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Table 3. Outflow and magnetic field alignment. The columns are (1) the region name, (2) the variance in polarisation percentage, (3) 
the average polarisation percentage measured across the region, (4) the mean position angle of the magnetic field vectors, (5) the outfiow 
direction, (6) the difference between the mean position angle of the magnetic field vectors and the outflow, (7) the opening angle of the 
outflow, (8) the inclination (from the plane of the sky) of the outflow, (9) notes and reference. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Region 


P 




P.A. Outflow P.A. 


m 


Opening Angle 


Inclination 




% 




O 


O 


O 






Cepheus A 


1.7 


1.3 


52 


45 


7 


15-20° 


pos (62°) 










90 


38 


60° 


pos (62°) 


DR21(OH) 


1.3 


1.0 


98 


90 


8 




pos 


GGD27 


1.2 


0.5 


105 


20 


85 


1° 


~pos 


GL2136 


2.3 


1.3 


102 


135 


33 


60° 


pos 


GL2591 


1.0 


0.6 


96 


45 


51 


<90° 


45° 










95 


1 


<90° 


45° 


GL437 


1.0 


0.6 


50 





50 


low 


pos 


IRAS 20126+4104 


0.8 


0.4 


74 


117 


43 


21° 


pos 










171 


83 


70° 


pos 


IRAS 20188+3928 


0.9 


0.3 


90 





90 




los 










45 


45 




los 


L1287 


1.4 


0.02 


75 


45 


30 


50° 


pos (60°) 


MonR2 


1.2 


0.7 


64 





64 




los 










45 


19 




pos 










90 


26 




los 


NGC6334A 


0.9 


0.8 


96 





84 




pos (10°) 


RCrA 


2.5 


7.0 


95 


90 


5 


60° 


pos 


S140 


1.5 


1.0 


70 


20 


50 




pos 










160 


90 




los 


S146 


1.3 


1.1 


86 





86 




los 


W49NW 


2.2 


3.2 


119 





61 


30° 


los 



(9) 

Notes/ References 



Patel et al. (2005) 
Patel et al. (20^) 
Lai et al. f 2003b) 



Gomez et al. f2003) : Martf et al. (1999) 
Kastner et al. (1994) 
jJut^warakori^_^_Cohen (2005 ) 
more coUimated than the 45° outflow 
Meakin. Hines fc Thompson (2005) 
Hofner et al. (2007) 
Hofner et al. (2007) 



Umemoto et al. (2000) 
Xu et al. (2006) 



Sarma et al. (2000) 
Anderson et al. (1997 ) 
Preibisch fc Smi th (2002) 



Ward-Thompson (1995 



Wu et al. (2005 



alignment of T-Tauri sta rs with the local magnetic field 
l|Menard fc Duchene|[20o3 ) reveal a possible connection be- 
tween the strength of the CTTS jets and their orientation 
with respect to the magnetic field. Interestingly, they con- 
clude that the CTTS's with jets align to the magnetic field, 
but as a whole sample (i.e. both CTTS's with and without 
jets), the population is randomly orientated with respect to 
the magnetic field, which they say suggests either the influ- 
ence of the magnetic field is dominant at large scales (whole 
cloud) but largely decreases on the much smaller scale of 
individual objects, or, the orientation of the CTTS's has 
changed since they first formed. 

Table [3] shows the mean polarisation percentage, posi- 
tion angles of the magnetic field and outfiow, and the differ- 
ence between the magnetic field and outfiow directions. The 
determined magnetic field vectors are not true (i.e. undirec- 
tional) vectors. They have a 180° ambiguity, and as such 
have position angles of between 0° and 180° . The magnetic 
field direction for each region is determined by calculating 
the weighted mean of the measured vectors. The smallest 
difference between the magnetic field direction and the out- 
flow direction is assumed, and the results are plotted in fig. ID 
Figure [2] is a cumulative distribution function, which shows 
that, for the whole sample (excluding S157, which as yet 
has no identified outflow, to the authors knowledge), given 
the (weighted) mean position angle of the magnetic field 
vectors, the magnetic fields within the sample appear ran- 
domly oriented with respect to the the jet/outflow direction. 
The Kolmogorov-Smirnov test reveals the whole sample has 
a 84.9% chance of being randomly orientated. 

The sample has been split into two sub-samples - those 



which have the outflow axis mainly in the plane of sky (i.e. 
i > 45° to the line of sight), and those which have the out- 
flow axis close to the line of sight (i < 45° ). GL2591 is not 
included in either of these sub-samples however, as the out- 
flow axis lies close to i ~ 45° . Whilst all of the sources with 
very small differences are in the plane of sky subsample, 
the cumulative distribution function does not deviate sig- 
nificantly from a random distribution, and the Kolmogorov- 
Smirnov test leads to a 57.2% chance of random orientation. 
However, the line of sight sub-sample, with no sources show- 
ing good alignment between magnetic fleld and outflow axis, 
does deviate significantly from the random distribution (a 
36.6% chance of random orientation from the Kolmogorov- 
Smirnov test), suggesting a relation between the misalign- 
ment of the (plane of the sky) field and outflow when the 
outflow lies close to the line of sight. If an alignment be- 
tween the magnetic fleld and the outflow axis does exist and 
the outflow is close to the line of sight, it would be difficult 
to assign a direction to the magnetic fleld, and so misalign- 
ment between the inferred fleld and outflow direction is more 
likely. 

Previous mid-infrared spectropolarimetry (|Aitken et al.l 

Il993l ) reveal that for their sample, overall, there is a large 
toroidal magnetic field component within the molecular 
structures associated with embedded young stellar objects. 
Their sample includes three sources from the sample in this 
paper - S140, MonR2 and GL2591. Their resuhs for S140 
and MonR2 agree with the findings of this paper - namely 
that there is a ~ 40 ° difference in alignment. However, they 
flnd the magnetic field of GL2591 to be perpendicular to the 
90° outflow, whereas we find it in almost perfect alignment. 
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This may be due to either a changing magnetic field mor- 
phology on diflFerent scales, or the magnetic field may have a 
stronger toroidal component closer in to the protostar/disc. 
Compared to the submillimetre, the mid-infrared emission 
comes from closer in towards the source, and it is higher 
resolution (has a smaller beamsize). 

The alignment analysis carried out here needs to be 
treated with caution, for some regions, e. g. MonR2, S140, 
the mean position angle may not accurately represent the 
field direction. MonR2, as with some of the other sources, is 
an extended source, and the polarisation pattern indicates 
abrupt changes in the direction of the magnetic field. For 
sources such as these, the best alignment comparison may be 
made using modal position angle values (in sufficient bins), 
or possibly from the map itself. 

Included in tableO where possible, are the opening an- 
gles for the jets/outflows. Upon quick inspection, it is clear 
that, with the information present, there appears to be no 
relation between the degree of coUimation and the align- 
ment. GL2591 has a well coUimated jet (position angle ~ 
95°) which is almost perfectly aligned to the magnetic field, 
whereas RCrA has an opening angle of 60° , and again is 
well aligned to the magnetic field. GGD27 has a well coUi- 
mated jet with an opening angle of only 1° , but the axis is 
almost perpendicular to the magnetic field. More informa- 
tion regarding the opening angles of such jets and outflows 
is needed to enable a more detailed analysis. 

5.2 Field Morphology & Submillimetre Emission 

A variety of continuum emission morphologies have been 
observed within this sample. Whilst the majority (11 out of 
16) of the sources have a seemingly roughly spherical dense 
core (marginally resolved in the JCMT beam), several of 
these sources also exhibit more extended emission associ- 
ated with the dense core (e. g. Cepheus A, DR21(0H), IRAS 
20188+3928). Of these sources, the polarisation percentage 
often decreases with increasing intensity across the cores, 
whereas in the extended emission the polarisation percent- 
age often remains stable. This may be due to the gas and 
dust being less dense in the regions of the extended emis- 
sion, or may be the magnetic field has a simpler, less twisted 
morphology in these areas. Morphologically, field lines both 
parallel (Cepheus A, W49) and perpendicular (DR21(OH), 
S157) to the major axis of the extended emission are ob- 
served. 

This sample also included dense cores which deviate sig- 
nificantly from spherical (at least in the plane of the sky, e.g. 
GL437, MonR2, NGC6334A). Of these sources the magnetic 
field is observed to be perpendicular to the major axis of the 
core in GL437 and RCrA. The field morphology of MonR2 
and NGC6334A are complex, with the fields changing direc- 
tion abruptly. This may be due to multiple cores unresolved 
by the JCMT beam. 

5.3 Field Strength & Cloud Support 

The plane of the sky magnetic field strengths have been 
calculated based on signal-to-noise (3) clipped vectors. The 
fields range from <0.1 mG (W49N) to 5.7 mC (Cepheus 
A), although the majority of them are around 0.2-0.4 mG. 



Out of the three strongest magnetic fields (Cepheus A, 
DR21(0H)N, RCrA), both DR21(OH)N and RCrA have 
field morphologies that may indicate ambipolar diffusion 
(the field is perpendicular to the major axis of the emission) . 
Cepheus A, however, has the strongest field, and the field 
is parallel to the major axis of the extended emission. Pre- 
vious polarimetric studies of W48 and S152 (jCurran et al.l 
l2004h . reveal magnetic field strengths and polarisation pat- 
terns consistent with the findings in this paper. W48 ex- 
hibits depolarisation across the main core, whereas the po- 
larisation percentage remains the same across the candidate 
HMPO W48W. S152 is a more complex region, consisting of 
ridges of gas and dust and the S152SE core - the candidate 
HMPO. The magnetic field is both parallel and perpendic- 
ular to the ridges in places, but across the S152SE core, the 
polarisation percentage remains at ~ 8% and the field is per- 
pendicular to the major axis of emission. The field strengths 
across the candidate HMPOs were calculated to be 0.7 mG 
for W48W and 0.2 mG for S152SE, consistent with the mag- 
netic field strengths calculated for the star forming regions 
in this paper. 



6 CONCLUSIONS 

We present the largest (to date) sample of high mass star 
forming regions observed using submillimetre polarimetry. 
We describe an improved method of SCUBA polarimetry 
data reduction. The magnetic field geometries (from the in- 
tensity weighted vectors, projected onto the plane of the sky) 
are presented, along with calculations of the core masses, 
densities and field strengths. 

The sources observed reveal a variety of morphologies 
for both the continuum emission and the magnetic field. 
The majority of the sources have dense cores that are unre- 
solved in the beam, but some also have extended emission 
associated with the cores. In several of the regions, the po- 
larimetry is uniform, suggesting ordered, relatively strong 
magnetic fields. A decrease in polarisation percentage across 
the cores is often seen, suggesting twisting or non-alignment 
of the dust grains in these cores. The polarisation percent- 
age often remains stable across the extended emission, which 
may imply that in these areas, the magnetic field has a sim- 
pler, less twisted morphology, or it may be that as the ex- 
tended emission is less dense, the dust grains do not become 
misaligned to the magnetic field. We see field morphologies 
perpendicular to the major axis of the continuum emission 
(DR21(0H), RCrA, S157), which could suggest ambipolar 
diffusion for these sources. It should also be noted that both 
of these regions which have outflows, have magnetic flelds 
that align well with the outflow axes. We also see fields that 
are parallel to the extended emission, e.g. Cepheus A, how- 
ever, the field in this region also aligns well with the outflow 
axis, and has the strongest magnetic field strength in this 
sample. It is interesting in this case that the extended emis- 
sion is in the direction of the outflow. 

An analysis of the mean position angle of the polari- 
sation vectors and the outflow axes has been carried out. 
Whilst the sample as a whole has a cumulative distribution 
function similar to that expected if the magnetic field and 
outfiow axes were randomly orientated, if the sample is bro- 
ken down by inclination, into those predominantly in the line 



Magnetic Fields in Massive Star Forming Regions 19 




^ — ' 1 1 1 1 1 1 1 1 1 

10 20 30 40 50 60 70 80 90 

Misalignement of magnetic field and jet/outflow axis [|d(P.A.)|] 

Figure 2. Cumulative distribution function of the difference is position angles between the weighted mean of the magnetic field vectors 
and the jet/outflow axis. The long dashed histogram is for all sources with outflow axes close to the plane of the sky (i > 45° ), the small 
dashed histogram is for all sources with outflow axes close to the line of sight {i < 45° ), the solid histogram is for the whole sample. 
The dotted line is the function expected for an infinite randomly oriented sample. 



of sight, and those predominantly in the plane of the sky, 
the Une of sight sample distinctly favours non-alignment, 
which would bo expected if the magnetic field had a large 
line of sight component too, as it is then difficult to establish 
a magnetic field direction from the polarimetry. The plane 
of the sky sample only shows a marginal increase in align- 
ment though. We discuss the caution needed in interpreta- 
tion of this analysis as the mean position angle of the vectors 
may not be the best representation of the field direction for 
all sources. In the sources where the field direction changes 
abruptly, modal averages (in reasonable bin sizes) may be 
more representative. Also, the polarimetry maps remain a 
very good way of analysing the alignment. 
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